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We study theoretically the quantum size effects in a one-dimensional semimetal by a Boltzmann 
transport equation. We derive analytic expressions for the electrical conductivity, Hall coefficient, 
magnetoresistance, and the thermoelectric power in a nanowire. The transport coefficients of 
semimetal oscillate as the size of the sample shrinks. Below a certain size the semimetal evolves 
into a semiconductor. The semimetal-semiconductor transition is discussed quantitatively. The 
results should make a theoretical ground for better understanding of transport phenomena in low- 
dimensional semimetals. They can also provide useful information while studying low-dimensional 
semiconductors in general. 
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I. INTRODUCTION 

Quantum size effects (QSEs) arise if the magnitude 
of at least one dimension of the sample is comparable 
to the de Broglie wavelength of carriers in the mate- 
rial. Experimental study of such effects started by an 
investigation of semimctallic bismuth thin films a long 
time ago [l|, followed by a theoretical model to explain 
the measured data soon after 0]. Since then the field 
has preserved its freshness to date, expanding its do- 
main to include the study of size effects in conventional 
metals Q, semiconductors 4], superconductors [jl and, 
since recently, nanoclusters 6[ , carbon nanotubes [7| and 
fullerenes Q. Recent advances in nanofabrication tech- 
niques have made it possible to extend experimental 
investigation of QSEs in various types of semimctallic 
structures from two-dimensional (2D) thin films to one- 
dimensional (ID) nanowires @-[I3]- On the other hand, 
despite the advances on an experimental front, there is 
not a theoretical model for quantitative understanding 
of measured data in such structures. Here, we general- 
ize the theoretical 2D model of Sandomirskii to study 
the confinement phenomena in a ID regime. Emphasis 
will be on QSEs in semimetals and on the semimetal- 
semiconductor (SM-SC) transition. To illustrate the ap- 
plicability of the model, size dependence of some of the 
transport coefficients in bismuth nanowires will be ad- 
dressed. 



II. SEMIMETALLIC REGIME 

In a semimetal the conduction and valence bands over- 
lap by a value A (see Fig. 1). At low temperatures, 
k^T <C A, as the sample size shrinks this overlap de- 



creases, finally reducing to the separation of the bands 
and the formation of an energy gap £ g . This is an imme- 
diate consequence of an alteration in the energy density 
of states. Let us consider a wire with dimensions w, t, 
and L (width, thickness, and length, respectively). The 
single-particle wave functions are given by [l8j 
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The corresponding energies for electrons are 
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where we have used the shorthand notation 
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Using the expression for electron density of states 
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and the Fermi-Dirac distribution for electrons with chem- 
ical potential n e , f e (£) = [1 + exp(£ — fi e ) / 'k^T]^ 1 , 
through the Sommerfeld expansion, we get for the vol- 
ume concentration of electrons at low temperatures 
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Above, O is the Heaviside step function and s is the spin 
degeneracy The derivation of an analogous expression 
for the concentration of holes p(£) is straightforward. 
The electron and hole chemical potentials fi e and /J 1 , 
measured from the bottom of a conduction band and the 
top of the valence band with respect to the Fermi level /i, 
can now be obtained from the charge neutrality condition 
n(£ ) = p(£). Writing the Fermi energy as fi e = fi% + (ly, 
we arrive at /x^/e^ = fi^/e 1 ]. Here j = x,y. Moreover, 
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FIG. 1: Color online, left) Overlapping bands in a semimetal. 
right) Dependence of the carrier concentration n, normalized 
by that for a wire with w/wo = 20 and t/to ~ 20, on the 
reduced width and thickness. A x — 0.7A. The bismuth wire 
is aligned along the bisectrix axes. The values for effective 
masses (in units of free-electron mass) are m% = 0.00139, 
m e y — 0.291 and m% = 0.0071 for electrons, and m£ = rriy = 
0.059, and m h z = 0.634 for the holes. 
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FIG. 2: Color online. Contribution of electrons to the elec- 
trical conductivity (triangles) and to the Seebeck coefficient 
(circles) of a bismuth wire versus its width at a fixed thickness 
t/to = 4.1. (This corresponds to a thickness of about to ~ 26 
mn). Other parameters are as given in Fig. 1. 



using the relation A = /i e + fi h , we get for the partial 
chemical potentials 



J = x,y. 
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Above, the energy overlap was represented as A = 
A x + A y . (Physically, A x and fi% can be interpreted 
as x components of the bands overlap energy A and 
the Fermi level fi e in the reciprocal lattice space; A^ 
is the overlap energy due to the confinement of carri- 
ers in the x direction.) As the sample dimensions be- 
come smaller, the conduction and valence bands slide up- 
ward and downwards respectively and the band overlap 
A gradually diminishes. The chemical potential /i, how- 
ever, remains intact. This is illustrated in Fig. 1. Let 
us assume that the band overlap vanishes at a width wo 
and at a thickness to, namely, e x (wo) + e x (wo) = A^ and 
e y(*o) + s y(to) = Ay. Using these relations together with 
Eq. (6), we get for the width and thickness at which the 
semimetal-semiconductor transition happens 
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Here Mj = m|m^/(m| + m!-) and j = x,y. The energy 
gap can now be written as £ g = A x (^-) 2 + A y (if) 2 - A. 
Designating r w = w/wq and r t = t/to, with 

= 1 - {i/r w ) 2 + {Hi/ nl) [1 - (j/r t f] , (8) 

we find for the normalized electron concentration 
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Above, [x] stands for the integer part of x. Figure 1 
depicts the dependence of reduced electron concentra- 
tion on width and thickness. Quantum size effects reveal 



themselves as steps in carrier concentration. By reduc- 
ing the sample size the electron density reduces in steps 
ultimately becoming negligible below transition width 
(thickness). Within each step, however, the electron den- 
sity varies nonmonotonously reaching its maximum value 
at a certain point. The positions of these maxima in each 
direction depend strongly on the effective masses, i.e., on 
the crystal structure of the wire. This is of vital impor- 
tance from an experimental point of view; observation of 
the effects discussed below presumes a well-defined crys- 
tal orientation and due care must be paid to the fabrica- 
tion of wires with comparable microscopic structures. 

In what follows, we will utilize the Boltzmann trans- 
port equation (BTE) to derive various kinetic coefficients 
of interest [l9| . Let us mark the unperturbed and per- 
turbed distribution functions with fo and /, the carrier 
charge with q £ {q e ,Qh}, and the velocity operator with 
hv(k) = Vk£(k). In the presence of an external electric 
field E, the BTE reads 



dtf + v ■ V r / + (g/h)E ■ V k / = (d t f) s 
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Through the linearization of this equation and the intro- 
duction of a relaxation time t, the kinetic coefficients of 
a ID system can now be obtained from 

(L n ) aS = ^ J(-dJo)r(k)(£(k) - ^)"v Q v 5 dk, (11) 

where Lo = <J stands for the electrical conductivity and 
S = L\/(qTLo) is the thermoelectric power (Seebeck co- 
efficient). To obtain an expression for the relaxation time 
we assume that TV scatterers, each with a strength Vb, 
are randomly distributed at positions R_, along the wire, 

^( r ) = Y^j=iVo& ( r — an d make use of the Fermi's 

golden rule = {2ir/h)\{i | V(r) \ f)\ 2 g(£ f )tt. Here 

= wtL is the sample volume. Averaging over the con- 
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figuration of scattering centers yields 
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volume density of scatterers. Electronic contribution to 
the electrical conductivity now reads 
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Here the outer summation accounts for the contribution 
of different subbands (m,n). Similarly, one obtains an 
analogous expression for the hole conductivity a h . The 
total electrical conductivity is the algebraic summation of 
electron and hole contributions. Figure 2 illustrates the 
contribution of electrons to the reduced electrical conduc- 
tivity CT e /a^ Q as a function of the wire width w/wq. The 
Seebeck coefficient can be evaluated using the Cutlcr- 
Mott relation valid for a degenerate gas at low temper- 
atures 0, S a5 = ^-^f^\n[g(£)v a (£)v s (£)T(£)}\ £ ^. 
Substituting for the corresponding quantities gives us 
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Dependence of the total electronic thermopower S e = 
E TO n Smn on the wire width is shown in Fig. 2. The 
thermopower observable in the measurements, covering 
the contribution of both the electrons and the holes, is 

cr"S"+cr h S h 



given by S 



III. SEMICONDUCTING REGIME 

In the semiconducting regime the carrier gas becomes 
ultimately nondegenerate and the distribution function 
will obey the Maxwell-Boltzmann (MB) statistics. De- 
noting the distance of the Fermi level to the bottom of the 
conduction band with [i e and to the top of valence band 
with n h , the distribution function can now be expressed 
as f j {£) w e -P(£ 3 +^) w ith (3 = l/k B T and j = e,h. 
Electron and hole concentrations for the intrinsic case, 
n = p (per unit volume) , are now obtained from 
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with the gap energy given by £ g = yf + fj, . The chemi- 
cal potential jj,, through the charge neutrality condition, 
reads 
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Generally, the MB distribution makes the derivation of 
analytic expressions for the kinetic coefficients more dif- 
ficult and in most cases one has to resort to numeri- 
cal methods [22j|. An alternative way, especially suit- 
able at lower temperatures, would be to estimate the 
quantities of interest with their average values. For the 
average value of the relaxation time over energy distri- 
bution, (t) = fT£g(£)d £ fd£ j f£g(£)d e fd£, assuming 

(l/vmn)) — Em'n' (l/( r mn— >m'n' 

)), we get 
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Substituting for the relaxation time in Eq. (11), the elec- 
trical conductivity follows: 
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The thermopower can be evaluated correspondingly, 
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At low temperatures, in the absence of intersubband scat- 
tering, one can limit himself to the contribution of the 
first subband only, that is, A™? n , = A\\ = 9, £f } w £{ x 
and E£j ~ £'i 1 . As a consequence, the electrical con- 
ductivity has exponential dependence on the gap energy 
£ g (and thus on the wire dimensions) while the ther- 
mopower does not. It would be in order to compare the 
expression above to that for a three-dimensional semicon- 
ductor with cubic symmetry given in Ref . fl9l] . namely, 
S e = ^-(f + c + (3fi e ). Here, c mirrors the energy depen- 
dence of the relaxation time, r e (£) oc £ c . 



IV. GALVANANOMAGNETIC COEFFICIENTS 

Let us assume that the wire is located in an external 
magnetic field B = (B x ,0,0) perpendicular to the wire 
axes and an external electric field E z is applied along 
the wire. The cyclotron frequency is marked with ui c . 
Moreover, let the cyclotron radius be much larger than 
the lateral dimensions of the sample, r c 3> w,t. The 
Hall coefficient R and the transverse magnetoresistance, 
= p< " B \ n ^ Q,> with p = 1/cr, can be easily evaluated by 
employing the two-carrier model [19|, |22j . We obtain for 
the Hall coefficient 
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((T e + o h ) 2 + (a e a h ) 2 (R e + R h ) 2 B 2 



(20) 



4 




w/w 



FIG. 3: Color online. Dependence of the Hall coefficient (tri- 
angles) and the magnetoresistance (circles) on the width in a 
weak magnetic field. Contributions from both the electrons 
and the holes, as given by the two-carrier model, are taken 
into account. 

and for the magnetoresistance 

_ Ap _ a e a h (u e +u h ) 2 B 2 

V= ~p~^~(<7 e + a' 1 ) 2 + (a e u h - a h u e ) 2 B 2 ' ' ' 

Above, i? J = (nqj)^ 1 with j = e,h is the single-particle 
Hall coefficient and vP stands for the carrier mobility. 
Figure 3 illustrates the dependence of the Hall coefficient 
and magnetoresistance in weak magnetic fields (co c t <C 
1) on the wire width at a fixed thickness. Below, we 
compare some of the findings of our study with the recent 
experimental data and with the other existing models. 



V. DISCUSSION 

First, let us compare our model to a recent theoretical 
study of the SM-SC transition in cylindrical nanowires 
in which the electron and hole energies are expressed in 
terms of the Mathieu (or Besscl) functions [l3[ . This 
has to be compared to the expressions in Eq. (7). While 
Ref. [HI predicts the enhancement of QSEs as a conse- 
quence of anisotropy, our model foresees sensitive varia- 
tion of the SM-SC transition width (and thickness) and 
alteration in the form of carrier density and transport co- 
efficients as the values of effective masses change (Figs. 
1-3). Furthermore, the model in Ref.[l3[ misses the is- 
sue of transport. As it comes to practical applications, 
in most cases, the (lithographically fabricated) wires are 
of rectangular cross section as modeled in our study. To 
make a quantitative comparison between predictions of 
our model and the other existing ones, we consider a bis- 
muth wire aligned along the bisectrix axis with effective 
masses given by Ref.[2fJ and presented in Fig. 1. The 
waveguide-analogous model in Ref. [l3[ predicts a critical 
diameter of d c « 32.5 nm for the SM-SC transition. An 
earlier similar model, instead, gives d c m 48.5 nm [l2T |. 



Correspondingly, in our model, substituting for the val- 
ues of effective masses in Eq. (7) and with A x = 0.7 A, we 
obtain wq ~ 102 nm and to w 26 nm. Now, the diameter 
of an equivalent cylindrical wire, i.e. a wire with an equal 
cross section, can be estimated as d c = \J 4iuo^o/ 7r ~ 58 
nm. 

As is the case with the electrical conductivity and the 
Seebeck coefficient, at a fixed thickness, all the oscilla- 
tions have the same periodicity given by Eq. (7). With 
chosen parameters given above, this yields wq ks 102 
nm. The corresponding value observed in the experi- 
ments with a 2D n-type bismuth film, given in Ref. fl4| . 
is about 100 nm. It is also a remarkable fact that, as pre- 
dicted by the theory, the measured value of film thick- 
ness for the semimctal-semiconductor transition equals 
the magnitude of periodicity near transition. However, 
in contrast to the predictions by our model, in experi- 
ments the period of oscillations becomes smaller as the 
transition thickness is approached. Also, the measured 
amplitudes of oscillations show less regularity than those 
envisaged by the model. While part of the discrepancies 
can be attributed to the structural nonidcalities in the 
grown films, the others can be due to uncertainties in 
the chosen values for the effective masses and due to the 
fact that the nondiagonal contributions in the effective 
mass tensor were ignored. Naturally, it would be more 
desirable to compare the theoretical predictions to the 
experimental data in a single one-dimensional wire with 
a well-defined crystal orientation. To our knowledge, de- 
spite recent experimental efforts with solitary nanowires 
of bismuth [TH, [l6[ , the oscillatory dependence of kinetic 
coefficients on wire dimensions are yet to be observed. 
For the verification of predicted effects, we consider it 
crucial to have nanowires with comparable morphology, 
and to characterize them with the precession four-probe 
measurements. 

Dependence of the SM-SC transition on temperature 
and on the dopant concentration was investigated in a 
recent study of nanowires made of alloys of bismuth [l2| . 
There, a scheme based partly on the measurement of 
the Seebeck coefficient, was introduced to differentiate 
scmimetallic nanowires from semiconducting ones. Cor- 
respondingly, in our model this can be performed by 
comparison of the form of the measured thermopower to 
the forms of Seebeck coefficients in SM and SC regimes, 
Eqs. (14) and (19), respectively. At low temperatures 
the latter reduces to S e w |&(§ + (3p e + 08^) and de- 
pends inversely on temperature, the former instead has 
a linear T dependence. Also, at higher temperatures, 
the thermopower in semiconducting regime saturates to 
an asymptotic value, S e ~ % — , comparable to (and 

slightly smaller than) the bulk value S e w |^(f + c). 
All these predictions agree with the experimentally ob- 
served behavior of the thermopower illustrated in Fig. 3 
of Ref. [ll[ ■ Our model does not account for the effect of 
dopant concentration on the transport properties of the 
alloys. However, in principle, there is not any restriction 
to extend it further to cover the doping effects too. 



5 



Finally, let us briefly discuss the relevance of our study 
to the similar phenomena in carbon nanotubes (CNTs). 
Depending on its crystal structure, a CNT can be metal- 
lic or semiconducting. Typically the energy gap of a 
semiconducting tube is small, ~ 0.5 eV. Any disturbance 
(bending of the tube for example) can alter the crystal 
structure and introduce a larger energy gap or reduce it; 
a metallic tube may turn semiconducting or vice versa. 
An analogy to the quantum size effects discussed above 
is apparent. Accordingly, we may anticipate the stepwise 
change of carbon nanotube conductance as a function of 
mechanical disturbance. The implementation of such an 
idea could be of interest in nanoelectromechanical ap- 
plications; in designing an ultrasensitive transducer, for 
instance. 



dimensional semimetals is introduced. The model can 
be extended to include the nonparabolic contributions to 
the dispersion relation and to cover doped semimetals 
and semiconductors at finite temperatures as well. Con- 
sidering the nonmonotonic oscillatory dependence of the 
transport coefficients on the sample size, the analytic re- 
sults presented here can be readily utilized to optimize 
various figures of merit in different applications. Similar 
approaches to the one introduced here can be adopted for 
the study of size effects in other low-dimensional systems, 
such as quantum dots or mesoscopic superconductors, for 
example. 



VI. CONCLUSIONS 



To summarize, a model for the study of confinement ef- 
fects and transition to the semiconducting regime in one- 
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